Casimir energy in multiply connected static hyperbolic Universes 
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We generalize a previously obtained result, for the case 
of a few other static hyperbolic universes with manifolds of 
nontrivial topology as spatial sections. 

PACS numbers: 98:80Cq, 98:80Hw 

I. INTRODUCTION 

As is well known, Einstein equations (EQ) restrict the 
local geometry of spatially homogeneous and isotropic 
spacetimes, to those of R^, S^, or H^. Recent observa- 
tional data indicate that the curvature of the universe 
is small, without ruling out the case of negative curva- 
ture. On the other hand, the EQ are insensitive to a 
global nontrivial topology of space, which can be a com- 
pact hyperbolic 3-space A4, which is isometric to a quo- 
tient space /r. Here F is a nontrivial discrete group of 
isometrics ( known as holonomy group), which acts freely 
and properly discontinuously on the covering space H^. 
Also important is the fact that F is isomorphic to the 
fundamental group tti{M), which is a group of homo- 
topy classes of maps of the circle into M 0. Since 
Tri{M) is nontrivial, A4 is multiply connected. 

The 3-space Ai may be represented by a fundamental 
polyhedron FP in H^, with an even number of faces, 
whose copies 'y{FP), 7 e F, fill up the entire H^. The 
faces of FP are pairwise identified by the basic elements, 
or generators, of F. The resulting manifold is a bundle 
with discrete fibers Tp over base points p in the funda- 
mental polyhedron. These fibers are the points of the 
quotient space. 

Among the first applications of the topology consider- 
ations, there was an attempt to explain multiple quasar 
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images §]. 

For recent reviews of topology in connection with cos- 
mology, see Q, and the article |^ for compactifica- 
tions of the 3-sphere. 

The first astrophysical limits on the topology of the 
universe were obtained for a 3-torus . Accordance with 
the homogeneity of the CMBR puts a lower limit on the 
size of the fundamental cell, about 3000 Mpc, which is a 
cube in the cases of and [Q. Later on, it was shown 
that this result is very sensitive to the type of the com- 
pactifications of the spatial sections. For a universe with 
spatial sections x R, the fundamental cell's size is 
about 1/10 of the horizon, and is compatible with the 
homogeneity of the cosmic microwave background radia- 
tion (CMBR) §. 

In compact universes, the pair separation histogram 
would present spikes for characteristic distances. At first 
it was thought that this technique, known as the crystal- 
lographic method, was able to totally determine the topol- 
ogy of the universe It turned out that the crystallo- 
graphic method only applied when the holonomy group 
contained at least a Clifford translation, i.e. a transla- 
tion which moves all the points by the same distance |l^] 
and ||ll| . Generalizations of the crystallographic method 
were proposed, for example in [ p^ . 

Also in compact universes the light front of the CMBR 
interacts with itself producing circles in its sky pattern 

ill. 

A recent result has called our attention to the possi- 
bility that methods based on multiple images will prove 
not to be efficient [Q. The reasoning is that, accord- 
ing to observations, the curvature is very small, so the 
fundamental regions are so big that there has not been 
time enough for the formation of ghost images. The re- 
sult is that for low curvature universes such as ours, only 
compact universes with the smallest volumes could be 
detected by pattern repetitions. 

A very attractive argument in favor of compact hyper- 
bolic manifolds is related to prc-inflationary homogeniza- 
tion through chaotic mixing [n5|. The effect is the same 
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that arises in compact hyperbolic surfaces. The geodesic 
motion on a surface of genus g > I shows the absence 
of KAM torus . Not only it is ergodic but also sat- 
isfies the Anosov property, which indicates the presence 
of strong chaos The chaotic properties of the g = 2 
torus were previously studied, for example in 



We extend our previously obtained result Il8| (see also 
[19 1) for a few more compact hyperbolic universes. In 
[18|, we numerically calculated the Casimir effect of a 
scalar field for a static universe whose spatial section is 
the Weeks manifold, the smallest volume (with curva- 
ture normalized to K = —1) compact hyperbolic mani- 
fold known. The outcome of our calculation is in fact a 
Casimir energy density pc- 

In compactifications of flat space, the strictly speaking 
Casimir energy can be analytically obtained through the 
very elegant zeta functions techniques; for a recent review 
see [Q. Unfortunately this beautiful formalism does 
not yield an analytical result in the case of hyperbolic 
compact manifolds, since the spectrum of the Laplace- 
Beltrami operator can only be determined numerically. 

We use the point splitting technique in the covering 
space, which is a static hyperbolic universe. The obtained 
propagator is exact and possesses information about the 
global properties of the manifold, in the sense that the 
infrared modes are taken into account. 

When the spacetime is multiply connected, the prop- 
agator is obtained as the usual sum over paths: all 
geodesies connecting the two points are taken into ac- 
count. 

We find a static hyperbolic solution for the EQ in sec- 
tion II. In section III, we write the expression for the vac- 
uum expectation value of the energy-momentum tensor 
in compact hyperbolic universes. We obtain in section 
IV the numerical values of the Casimir energy density in 
a few multiply connected static spacetimes. Our conclu- 
sions are presented in section V. (We use natural units, 
G = c = h = 1, except in section II.) 



II. QUANTUM FIELD THEORY IN THE STATIC 
UNIVERSE Rx 

The hyperbolic space sections H^, can be realized as a 
surface 

(x - x'f + {y- y'f + {z- z'f ~ {w - w'f = -a", 

(1) 

imbedded in a Minkowski 4-space 

df = dx^ + dy^ + dz^ - dw'^ . 

As this space is homogeneous, we explicitly write the ori- 
gin of coordinates (x', y',z\w'). It can easily be seen that 
its isometry group is the proper, orthocronous Lorentz 
group S'O^ (1,3), which is isomorphic to PSL{2,C) = 



SL(2,C)/{±I} I|T|. With the constraint of Eq. (|T|) on 
the line element we obtain 



dx^ + dy^ 



{{x — x') dx + {y — y') dy + (z ~ z') dzY 
{x - x'f + {y- y'f + {z- z'f + 
ds^ = -dt^ + df = g{x,x')f_,t,dx'^dx'', 



(2) 



where we interchangeably write (x*^, x^, x^, x^) < — s- 
(t,x,y,z). Both connections, V^; and V^;', compatible 
with the metric of Eq. (|^), can be defined through 



V^g(a;, a;')a/3 ^ 
Wf,>g{x,x')a0 



, 

E 



(3) 
(4) 



The expression in Eq. (^ is the popular Robertson- 
Walker line element, which written in the Lobatchevsky 
form reads 



ds^ = 
with 



-dr 



sinh X 



[dx^ + sinh^ X {dO^ + sin^ &d(l)^)] , (5) 



{x^x'f + {y-y'f + {z-z'f 



As is well known, the EQ for the homogeneous and 
isotropic space sections in Eq. (|^), with a = a{t), reduce 
to the Friedmann-Lemaitre equations 



SttG 



A 

a" 



TT = -SnGp + A , 



where the right-hand side comes from the classical 
energy-momentum source for the geometry, T'^'^ = (p + 
p)u^u^ -\-pg^^ , plus the cosmological constant term Ag^"^. 

We assume that the universe was radiation dominated 
near the Planck era, hence p = p/3, and obtain the fol- 
lowing static solution 




SttG 
ds^ = -dt^ 



[dx^ +sinh2x(d6|2 +sin2e'd</>2)] , (6) 



where the cosmological constant is negative. 

We now wish to evaluate the vacuum expectation value 
of the energy density for the case of a universe consisting 
of a classical radiation fluid, a cosmological constant, and 
a non-interacting quantum scalar field 4>. The solution of 
EQ is given in Eq. (|^) , where the quantum back reaction 
is disregarded. We use the point splitting method in the 
universal covering space R x , for which the propaga- 
tor is exact. The point splitting method was constructed 
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to obtain the renormalized (finite) expectation values of 
the quantum mechanical operators. It is based on the 
Schwinger formalism p2| , and was developed in the con- 
text of curved space by DeWitt |^ . Further details are 
contained in the articles by Christensen |^^, For a 
review, see p6| . 

Metric variations in the scalar action 



S 



where T is the time-ordering operator. Taking its real 
and imaginary parts, 



G{x,x') = Gs(x,x') + -G^^\x,x') 



(11) 



we get for the Hadamard function 

G'^^\x,x') = (O|{0(x),0(a;')}|O) = 2ImG(a;,a;') . 



with conformal coupling ^ — 1/6, give the classical 
energy-momentum tensor 
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GauC, 



-m^gp^cj)^ , (7) 



where G^i, is the Einstein tensor. As expected for mass- 
less fields, it can be verified that the trace of the above 
tensor is identically zero if m = 0. Variations with re- 
spect to (f) result in the curved space generalization of 
the Klein-Gordon equation. 



R 

□ 0-771^,^ = 

6 



(8) 



The renormalized energy-momentum tensor involves 
field products at the same spacetime point. Thus the 
idea is to calculate the averaged products at separate 
points, X and x' , taking the limit x' — > a; in the end. 



^\Tpv{x) |0) = lim T(x,x')p^ , 



(9) 



with 
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(V^V, + V,,. V.O + j^gix)^,. (□ + □') 
T2 (^(^^^'^ ~ \^(^)9{x)pi.j - ^m^g{x)p^ 



(10) 



G^^\x,x') , 

where the covariant derivatives are defined in Eqs. (|^) 
and (III), and G'-^-' is the Hadamard function, which is 
the expectation value of the anti-commutator of 4>{x) 
and 4>(x') (see below). We stress that the quantity 
T(x, x')p_i, only makes sense after the limit in Eq. |9| 
is taken. The geometric quantities such as the met- 
ric and the curvature are regarded as classical entities. 
g{x)pi, — g{x^' — Q)fiu is obtained from the line ele- 
ment in Eq. (g). 

The causal Green function, of Feynman propagator, is 
obtained as 

G{x,x') = i(O|T0(x)0(x')|O) , 



III. THE FEYNMAN PROPAGATOR AND THE 
CASIMIR ENERGY DENSITY m RxM 

Green functions, as any other function defined in the 
spatially compact spacetime i? x , must have the same 
periodicities of the manifold M itself. One way of im- 
posing this periodicity is by determining the spectrum of 
the Laplacian, which can only be done numerically. 

Another method imposes the periodicity by brute 
force, 

fM{x) = ^ filx) . 

The above expression is named the Poincare series, and 
when it converges, it defines functions /x on the mani- 
fold M. 

We define the operator 



F{x, x') = F{x)/y/^5{x ~ x') , 



(12) 



where F{x) = □ — R/Q — m? , and introduce an auxiliary 
evolution parameter s and a complete orthonormal set of 
states \x), such that 



G{x,x') = {x\G\x') , 
F{x,x')^{x\F\x') , 



G = I I e 
/o 



(13) 



This last equation implies that G = {F — iQ) hence 
the causal Green function becomes 



G{x,x')—i / ds{x\exp{—isF)\x') 



(14) 



and the matrix element (a:| exp(— jsi^)|a::') = (a::(s)|a;'(0)) 
satisfies a Schrodinger type equation, 

^^Jxis)\x'm =(□-§- m^) (^(s)l^'(O)) . 

Assuming that {x{s)\x' {0)) depends only on the total 
geodesic distance —(t — t'Y + a?x^ i with the spatial part 
a?"X^ derived from Eq. (p|), the above equation can be 
solved, and we get 
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{x{s)\x'm 



—iX exp{iTO S + — — a^x ]/4s} 



sinhx 



(47rs)2 



(15) 



Substituting this solution for the integrand in Eq. (14) 
gives for the Feynman propagator 



G{x,x') 



m X 



mJ {t - t'Y - a^x^ 



Stt sinh x 



^{t - t'f - a^x^ 



(16) 



where H\ ' is the Hankel function of the second kind and 



order one. 

The Hadamard function can be obtained from Eqs. 
(H) and (0), 



G^^\x,x') 



m X 



Ki ( mJ - {t - t'Y + a^x^ 



2Tr^ sinh x 



{t - t' f + a^x^ 



(17) 



where Ki is the modified Bessel function of the second 
kind and order one. The massless hmit m = can be 
immediately be checked: 



G^'Hx,x')m=0^ 



X 



1 



27r2 sinh (x) \~{t~ t' f + a^x^ 
Remembering that for a — > cx3 



sinh x — o. \/ {x — x'f + {y — y'Y + [z — z'f x j 

the well known Minkowski result is recovered for the mas- 
sive and massless cases, 



G^'-\x,x') 



K, mJ~{t-t'Y + 



27r2 



G^'\x,x'U=o 



- (< - t'Y + r2 
1 



27r2 \-{t-t') 



where r is the geodesic distance in the spatial euclidean 
section. 

Substituting Eq. (1^) and the covariant derivatives (||) 
and (^) into Eq. (|lO), we obtain T{x,x')^^. 

The Klein- Gordon equation remains unchanged under 
isometries, 



□ -|-m2)0 



R 
6" 



in the Green functions over the discrete elements of the 
group r is well defined. 

In A4 — H^/T, a summation over the infinite num- 
ber of geodesies connecting x and x' is obtained by the 
action of the elements 7 G F, which are the generators 
gi and their products (except the identity; see below), 
on x' . Since F is isomorphic to TTi{Ai) - see section I - 
each geodesic linking x and x' in A4 is lifted to a unique 
geodesic linking x and ^x' in H^. Thus from Eq. (|) we 
get 

Pc^{0\T{x),,,\0)j^u^u-' 

= uV hm Y.^{x,^x') . (18) 

The infinite summation occurs because the spacetime 
R X A4 is static, so there has been enough time for the 
quantum interaction of the scalar field with the geometry 
to travel any distance. Since we know the universe is ex- 
panding, the infinite summation is not physically valid. 
The presence of the mass term, however, naturally intro- 
duces a cutoff. 

In Eq. ( p^ the subscript 7^1 means that the direct 
path is not to be taken into account. We shall show, 
following p6[ , that this exclusion is indeed equivalent to 
a renormalization of the cosmological constant. 

First remember that the effective action W is given by 



jw 



where the action for the scalar field can be transformed 
into a gaussian type after integration by parts, 

S = -\j d^x^/^) J d'y./^)4>{x)F{x,y)4>{y) . 

In a very informal way, the functional integration can be 
regarded as a usual gaussian integrald which, making use 
of Eqs. ( |T^ ) and (|l3|), results in 

.-1-1/2 
iF 



iW 



det — 
27r 



W = - ln[dct F] + const. 

SW = UrSilnF) = Ut[F-^SF] 



SW = -tr 
2 



SW 



W = — tr 
2 



e-''^SFds 





r°° p-i'^p \ 


tr 




. 2 \ 





ds ] ~ d^xC^tf 



where £^ is the Lie derivative with respect to the Killing 
vector ^ that generates the isometry, hence summations 



^This result is made more rigourous by assuming a particular 
complete representation basis for the operator F. 
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The trace of the operator is over the Hilbert space, de- 
fined in (|I3|), tvA = livax'-ixj ^/^d'^x{x\A\x'). Using 
the Schrodinger kernel {x\exp{—isF)\x'), given in Eq. 
(ITsI), we obtain the effective lagrangian 



(47r)2s3 



{x\ exp{—isF)\x') 



ds ^ AooV-g , 



ds 



(19) 



which shows that the direct path (7 = 1), corresponds to 
a divergent cosmological term Aqo- 



IV. CASIMIR DENSITY pc IN A FEW 
UNIVERSES 

According to quantum cosmology, a smaller universe 
has a greater probability of being spontaneously created. 
Also, the chaotic mixing becomes more significant for 
smaller volumes |l5| . We describe some spatially com- 
pact unive rses, with increasing volumes, in subsections 



[V A| - |[V D| . As seen in section I, manifolds M = H^/T, 
where F is a discrete subgroup of isometrics and is 
its universal covering, are multiply connected. 

The values of pc shown for each manifold were taken 
at points {6, ip) on the surface of a sphere inside its funda- 
mental region. For all of them the radius of the sphere is 
the same, d = ax = 0.390035. ..a, where d is the geodesic 
distance. Our result is displayed in FIGS. ^ ^, ^ and ^ for 
a scalar field with mass m — 0.4, and a metric scale fac- 
tor a — 10. Angles 6 and (p correspond to the co-latitude 
and longitude, so the lines 9 — and 6* = tt correspond 
to the poles of the chosen sphere. For each manifold we 
also write explicitly the radius of the inscribed sphere 
R 



inradius ■ 



Th e des cription that follows applies to all subsections 
' A - [ V D . The gi matrices that generate F were ob- 
tained with the computer program SnapPea [ pTf . The 
numerical code has been improved since our previous pa- 
per |18| . To yield a numerical result, the infinite summa- 
tion Wq) has to be truncated. Recall that this summation 
occurs in the covering space H^. We halted the summa- 
tion each time the action of the generators gi and their 



products on the origin (x^. 



^) = (0,0,0), reached a 



geodesic distance bigger than d = ax ~ 5. 29834. ..a. Care 
was taken so that no point was summed more than once. 
In other words, the summation (|l8|), which yelds pc, was 
truncated when the interior of the hyperbolic sphere of 
geodesic radius d = ax = 5. 29834. ..a, was covered with 
replicas of the fundamental region. We checked that ad- 
ditional contributions were unimportant for the evalua- 
tion of Pc- 



A. Weeks Manifold 

This manifold was discovered independently by Weeks 
psf and Matveev-Fomenko |^^, and is the manifold 
with the smallest volume (in units of a^) known, V = 
0. 942707. ..a^. Its fundamental region is an 18-face poly- 
hedron, shown in FIG. 0. The radius of the inscribed 
sphere is Rmradius = 0.519162. ..a. 

The vacuum expectation value of the 00-componcnt of 
the energy-momentum tensor, pc — T^^u^u'^ ^ as seen by 
a comoving observer, is shown in FIG. 0. 




FIG. 1. Fundamental region for Weeks manifold. 




FIG. 2. PC for Weeks universe. 



B. Thurston Manifold 



It was discovered by the field medalist William 
Thurston This manifold possesses a fundamental 

region of 16 faces, its volume \sV = 0.981369.. .a^ (FIG. 
I), and Rinradius = 0.535437. ..a. 

FIG. ^ shows the value of pc — T^yU^^u" — Too, as 
seen by a comoving observer. 
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FIG. 3. Fundamental region for Thurston manifold. 




On 271 



FIG. 4. PC for Thurston universe. 

C. Best Manifold 

This manifold was discovered as a by-product of a 
study of finite subgroups of 50(1,3) by a geometrical 
aproach pl| . Its fundamental region is an icosahedron 
with V = 4.686034.. .a^ and R^nrad^us = 0.868298.. .a, 
shown in FIG. |. 

The vacuum expectation value of pc — T^^u'^u'^ , as 
seen by a comoving observer, is shown in FIG. ||. 




FIG. 5. Fundamental region for Best manifold. 




It 2X 



FIG. 6. PC for Best universe. 

D. Seifert- Weber Manifold 

For this manifold, which was discovered by Weber and 
Seifert @, V" = 11. 199065.. .a^, R.^radtus = 0.996384.. .a, 
and the fimdamental region is a dodecahedron (FIG. |^). 

FIG. ^ shows the value of pc — T^j^^u^u'^, as seen by a 
comoving observer. 




FIG. 7. Fundamental region for Seifert- Weber manifold. 
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FIG. 8. PC for Seifert- Weber universe. 



V. CONCLUSIONS 

We explicitly evaluated the distribution of the vacuum 
energy density of a conformally coupled massive scalar 
field, for static universes with compact spatial sections 
of negative curvature and increasing volume: Weeks, 
Thurston, Best, and Seifert- Weber manifolds. As a spe- 
cific example, we chose m = 0.4 for the mass of the 
scalar field, and a = 10 for the radius of curvature. The 
values of the Casimir energy density pc on a sphere of 
proper (geodesic) radius d = 3.90035.. inside the funda- 
mental polyhedron for each of these manifolds are shown 
in FIGS, |[ ^, and |[ In all these cases it can be seen 
that there is a spontaneous generation of low multipolar 
components. As expected, the effect becomes weaker for 
increasing volume universes. 
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